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ABSTRACT 

We report on the detection of a hot subdwarf component in the Be binary 
system, 59 Cygni. The spectral signature is found in cross-correlation functions 
of photospheric model spectra with far-ultraviolet spectra obtained by the In- 
ternational Ultraviolet Explorer Satellite, and we used radial velocities from the 
cross-correlation functions to determine a double-lined spectroscopic orbit. The 
individual spectra of the binary components were extracted using a Doppler 
tomography algorithm. The flux of the system is dominated by the rapidly ro- 
tating Be star. However, the subdwarf contributes approximately 4% of the UV 
flux, and its spectrum bears a strong resemblance to that of the hot sdO star 
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BD+75°325. Based upon the appearance of the UV spectrum and the orbital 
elements, we present estimates for the stellar masses, radii, and temperatures. 
The presence of the hot companion causes excess emission from the outer part 
of the Be disk facing the companion. We present a set of red spectra that show 
the orbital phase variations of the He I A6678 emission formed in the heated re- 
gion of the disk, which probably occurs near the disk outer boundary. 59 Cygni, 
FY Canis Majoris, and <fi Persei comprise the known set of Be binaries with de- 
tected hot evolved companions, which are the stripped down remains of mass 
transfer. Their properties demonstrate that some fraction of Be stars were spun 
up through angular momentum transfer by Roche lobe overflow. 

Subject headings: stars: emission-line, Be — stars: individual (HR8047, HD 200120, 
59 Cyg) — stars: individual (BD+75°325) — stars: binaries: spectroscopic - 
stars: evolution — stars: subdwarfs 



Introduction 



The origin of the rapid rotat ion that characterizes the emission line Be stars is still 
unknown ( IPorter &: Riviniusll2003l ). One important clue is that Be stars are fr equently th e 
primary components of massive X-ray binaries with neutron star companions ( Reig boil ) 



and t hus, many are members of binaries that must have interacted in the past (IPols et al. 



f 99ll ). Past mass transfer from the neutron star progenito r to the mas s gainer probably 
resulted in a spin up of the gainer to near critical rotation (IPacketl Il98lh . and the current 
mass loss into the circumstellar disk of the Be star represents processes that shed some of the 
star's extreme angular momentum. In many cases, the donor star may lose enough mass to 
drop below the Chandrasekhar limit, and the result may be a hot helium star (spectral type 
sdO) or (eventually) a white dwarf remnant ( lGiesll2000f ). Such faint companions are very 
difficult to detect, but because they are hot, searches in the ultraviolet part of the spectrum 
are favored, and w e now have found the UV spectr al signatures of ho t helium cores in two 
Be binaries, <\> Per (jGies et al.lll998h and FY CMa JPeters et al-lboosh . 



Here we turn our attention to the remarkable collection of UV spectra in the archive 
of the International Ultraviolet Explorer (IUE) of the Be binary 59 Cygni (HD 200120, 
HR 8047, V832 Cyg). The star is the brightest componen t of a multiple syste m with a nearby 
companion Ab detected through speckle interferometry (IMason et al.ll2009l ). An even closer 
binary companion of the cent ral star Aa was suspec t ed for a lo ng time, and the f irst o rbital 
determinations were made by lTarasov fc Tuominenl (119871) andlRivinius fc Stefll ( 120001 ) who 
correctly identified the 28 d orbital period of the binary. lHarmanec et al.l ( 120021 ) presented 
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a joint photometric and spectroscopic analysis, and they determined orbital elements from 
radial velocity measurement s of t he Ha emissi o n line wings (formed in the disk close to 
the Be star). iMaintzl ( 120031 ) and iMaintz et al.l ( 120051 ) presented additional radial velocity 
measurements for the Be star, and they showed that anti-phase Doppler shifts are present in 
the He II A4686 absorption line, a fe ature only found in stars much hotter than the Be star 
primary (classified as B1.5 Vnne by lLeshlll968l ). Maintz and collaborators argued that the 
He II line forms in the atmosphere of a hot, faint companion and that this companion must 
be the stripped-down remains of the former mass donor star. They also documented how the 
hot companion illuminates the nearby rim of the Be star's disk and causes orbital-related 
emission line variations similar to those observed in the spectra of <j) Per and FY CMa. 

The spectral lines of a hot companion should be present in the ultraviolet spectrum, but 
they may be difficult to discern in individual spectra because of the complex line blending in 
the UV spectra of hot stars and the relative faintness of the companion compared to the Be 
star. The general appearance and long-term variations of t he str ong wind lines in the UV 
spectrum of 59 Cyg are described in detail by lDoazan et al.l ( 119891 ). Here we apply the cross- 
correlation method to dete c t the companion s pectrum and to measure its radial velocity 
variations (IGies et al.l Il998l ; iPeters et al.1120081 ). The measurements and the derivation of 
orbital elements are described in Section 2, and we use the radial velocity curves and a 
Doppler tomography algorithm to reconstruct the spectra of both components in Section 3. 
We describe in Section 4 how the flux of the subdwarf companion heats the nearby region 
of the Be star disk and creates orbital phase related emission line variations. We summarize 
these new results on 59 Cyg and their relation to other Be + sdO binaries in Section 5. 



2. Radial Velocities and Orbital Elements 

The archival spectra of 59 Cyg from observations made with the International Ultraviolet 
Explorer provide a valuable record of the orbital Doppler shifts associated with the binary 
motion. Here we describe the observational properties of the spectra, and then we outline 
our use of the cross-correlation method to measure radial velocities of first the hot secondary 
star and then the relatively cooler Be star. Because the radial velocity semiamplitude is 
much larger for the hotter star, we begin the fit of the orbital elements with the hot star 
velocities and then make a restricted fit of the associated elements for the Be star. We 
conclude this section with a comparison to earlier published work. 

We obtained the UV spectra of 59 Cyg from the NASA Mikulski Archive for Space 
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Telescopeg^. The selection focused on observations made with the IUE HIRES, Short Wave- 
length Prime (SWP) camera, which recorded the far-UV spectrum with a resolving power 
of R = 10000. There are a total of 194 such observations that were made between 1978 and 
1994. We further limited this sample to the 157 spectra with the best S/N ratio (usually 
> 10 in the better exposed parts of the spectrum). The echelle orders were extracted and 
combined using pipeline IUERDAF software, and then each spectrum was normalized to 
unity in relatively line-free regions and transformed to a standard, heliocentric wavelength 
grid in log A that spans from 1200 to 1900 A. The wavelength calibration of each spectrum 
was checked for possible offsets due to the star's position in the large aperture by compar- 
ing the positions of the strong interstellar lines with those in the global average spectrum. 
For each spectrum, we made cross- correlation measurements of the pixel shifts between the 
strong interstellar lines in a given spectrum with those in the average of all the spectra, and 
then we applied the median of these measurements to shift the observed interstellar spec- 
trum into alignment with the average interstellar spectrum. Then the interstellar lines were 
removed from each spectrum by replacement with a linear interpolation to the spectrum at 
the boundaries of each interstellar feature. 

Our primary goal is to detect the spectrum of the hot companion in order to measure the 
star's radial velocity. It is very difficult to measure the position of any particular line of the 
companion's spectrum because of the limited S /N of the observations and the small fraction 
of the companion's flux (see below). However, we can measure the signal of the ensemble 
of spectral lines of the companion by calculating the cross-correlation function (ccf) of the 
spectrum with that of a model for the hot star. This approach is successful because the 
companion is so much hotter than the Be star that almost all of its spectral features are 
comprised of higher ionization species, and consequently the ccf signal is only sensitive to the 
lines of the hot companion and not to those of t he cooler Be star. We co nstructed a hot star 



spectrum template from the spectral models of lLanz fc Hubenyl (120031 ) for T e g = 55000 K, 
logg = 4.4, Vsini = km s _1 , and solar abundances. The model was transformed to 
the observed wavelength scale by flux integration into the same log A wavelength system, 
and then the model was convolved with the instrumental broadening function. The ccf 
was calculated over the entire spectrum after first setting to unity the low wavelength end 
(1200 - 1226 A to avoid Lya) and those regions in the immediate vicinity of strong wind 
lines or interstellar lines. Consequently, the results are largely free from the influence of 
non-photospheric spectral features. 

We show the results of the ccf calculations in Figure 1 as a function of radial velocity 



1 http://archive. stsci.edu/iue/search.php 
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and orbital phase. The upper panel illustrates the difference between each ccf and the mean 
of all the ccfs in order to remove the broad and constant background ccf structure caused by 
low frequency correlation with the Be star spectrum. It is hard to find the peak signal from 
the spectrum of the hot companion in individual ccfs because the peaks have an amplitude 
not much larger than the noise. However, the signal is much more evident in the lower 
panel where the ccf differences are plotted as a gray scale intensity as a function of radial 
velocity and orbital phase. The bright, "S" -shaped curve reveals the orbital motion of the hot 
companion, despite the weakness of the signal, thanks to the large number of observations 
available. 

We estimated the radial velocities of the hot companion by fitting a parabola to the 
central ±20 km s -1 around the peak of each ccf. This was not always successful because the 
peak was too weak in many cases, and in the end we restricted the results to 132 spectra 
where a reliable measurement was possible. The results are summarized in Table 1 (given 
in full in the electronic version) that lists the heliocentric Julian date of mid-observation, 
the corresponding UT date, the IUE SWP number assigned to the spectrum, and orbital 
phase (see below). Columns 5-7 and 8-10 then give the measured radial velocity (V), 
uncertainty (a), and observed minus calculated residual (O — C) for both the primary Be 
star (see below) and hot companion, respe ctively. Note that in most cases that the formal 



uncertainty a (based in the formulation of IZuckerl 120031 ) is smaller than the absolute value 



of the residual O — C because the former does not account for uncertainties associated with 
the wavelength calibration of each spectrum. 

We then fit t he radial velocities to o btain orbital elements using the non-linear, least- 



squares solver of iMorbey fc Brosterhud (11974 ). The solution was obtained by assigning 



equal weights to each measurement, which is appropriate for the uniform S/N properties 
of the final set of measurements. The resulting elements and their uncertainties are given 
in column 4 of Table 2, where the subscript 2 is used to identify parameters for the hot 
companion. These include the orbital period P, the epoch of periastron T, the epoch of 
Be star superior conjunction T$c (equal to the epoch of inferior conjunction for the hot 
companion), the velocity semiamplitude K, the systemic velocity 7, the eccentricity e, the 
longitude of periastron oj, and the root mean square of the residuals from the fit. We adopt 
the standard orbital phase convention in which orbital phase zero corresponds to the time of 
periastron. The orbital velocity curve and measured radial velocities are plotted together in 
Figure 2. Note that the fitting scheme uses a version of the Levenberg-Marquardt method 
that estimates the parameter errors according to the covariance matrix and the final \ 2 °f 
the fit. 



We also used the cross-correlation method to measure radial velocities for the Be star 
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primary, but in this case we used the average of all the spectra as the spectrum template. 
This is a reasonable choice because the flux of the Be star dominates the far UV spectrum. 
The resulting ccfs show a very broad peak that reflects the large rotational broadening of 
the primary star's spectrum plus a weak and narrow central peak. We suspect that the 
latter component may result from residual interstellar lines, edges corresponding to cuts of 
the stronger interstella r features, and/or lines of the relatively stationary third component, 
Ab (jMason et al.ll2009l ). Thus, instead of fitting a parabola to the central part, we measured 
the ccf positions us ing the bisector of the wings that was calculated using the method of 
Shaffer et al.l ( 119861 ). The transformation from relative to absolute radial velocity requires an 
estimate for the velocity of the average spectrum. We found that the velocity offset between 
the average spectru m and a model (for T c s = 21750 K, logg = 3.8, Vsini = 379 km s _1 ; 
Fremat et aDhoOsh was 7.8 ± 7.4 km s 1 (by the wing bisector method). However, this ccf 
with the model was visibly asymmetrical, and the resulting bisector velocity depended on 
the part of the wing sampled. Consequently, the transformation from relative to absolute 
velocity remains problematical, so no correction was applied to the velocities of the primary. 
The resulting relative velocities, uncertainties, and residuals are listed in columns 5 through 
7 of Table 1 for all 157 spectra selected. 

A first fit of the primary velocities with only the period fixed resulted in estimates of 
e, oj\ = 0J2 + 180°, and T consistent within uncertainties with the respective values for the 
secondary but with much larger values of uncertainty. Thus, we simply fixed P, T, e, and u 
to those values derived from the larger amplitude radial velocity curve of the secondary, and 
then solved for K\ and 71 only. The derived parameters are listed in column 4 of Table 2, 
and the radial velocity curve and measurements are over-plotted with those for the secondary 
in Figure 2. With estimates in hand for the parameters of both components, we can then 
determine the mass ratio, Msin 3 i products, and the combined projected semimajor axis 
asinz, which are all listed in Table 2. 

Our results for the radial v eloci ty curve of the hot c ompanion are in excellent agreement 
with those from iMaintzl (120031 ) and iMaintz et al.l (120051 ). and this confirms their conclusion 
that the Doppler shifts of the He II A4686 line trace the motion of the hot star. On the other 
hand, our estimate for the semiamplitude of the Be star is much sm aller than that found by 
Maintz and collaborators, but agrees with the determination from lHarmanec et al.l (120021 ) 
based on the Doppler shifts of the Ha emission line wings. We think that t he diff erence arises 
because of residual emission in the He I lines measured by lMaintz et al.l (120051 ) . They point 
out that the emission component formed in the outer disk facing the hot companion (which 
is seen dramatically in He I A6678; Section 4) will probably partially fill in the absorption 
profiles of the other He I lines in such a way as to increase the apparent Doppler shifts. The 
UV lines, on the other hand, correspond to larger energy transitions than those of He I, and 



-7- 



the temperatures in the heated portion of the disk are probably insufficient to create any 
significant emission line flux in the FUV (and we see no evidence of such emission in the IUE 
spectra). The fact that the velocity semiamplitude K\ we derive from the UV photospheric 
spectrum is so similar to that from the Ha emission wings (formed in the inner part of the Be 
star's disk) is consistent with the assumption that the Doppler shifts we measure correspond 
to those of the Be star itself. Estimates of t he other orbital elemen t s (with t he po ssible 
exception of o j) are consistent with those from lHarmanec et al.l (120021 ) , iMainta (120031 ) , and 



Maintz et all (I2005[ ). 



3. Tomographic Reconstruction of the UV Spectra 



The nature of the faint companion star can be analyzed through an examination of its 
UV spectrum. However, because the companion is so faint compared to the Be primary, 
its spectrum cannot be easily ext racted from any indi vidual spectrum. Instead, we used 
a Doppler tomography algorithm (IBagnuolo et al.lll994l ) to extract the component spectra, 
and because this approach uses all the spectra, the reconstructed spectra have a S/N ratio 
improved by ~ y/n, where n is the number of spectra. We used the orbital solutions to 
determine the radial velocities of the primary and secondary for each spectrum, and we 
adopted a small flux fraction to begin (that was adjusted later to match spectral depths in the 
reconstructed spectra; see below). We adopted the same model spectra used for ccf analysis 
as starting values and we ignored the flux contribution from the faint Ab speckle companion. 
We ran the Doppler tomography for 50 iterations with a gain of 0.9, and tests showed that the 
resulting spectra were insensitive to all these assumptions. The reconstructed spectrum of the 
hot secondary shows considerable wander in the continuum on scales of tens of Angstroms, 
so we reset the continuum by forming Gaussian smoothed (FWHM = 610 km s _1 ) versions of 
both the reconstructed and model spectra and then by dividing the reconstructed spectrum 
by the ratio of the smoothed spectra. This method removes all long wavelength span features 
but retains the narrow features that characterize the spectrum of the secondary. The rectified 
version of the reconstructed secondary spectrum is illustrated in Figures 3 to 7 that show 
spectral subregions containing lines of the principle ions that characterize the spectrum. 



We compared the reconstructed spec tra to TLUSTY/ S YNSP EC models from lLanz fc Hubeny 
( 20071 ) for the Be star primary and from lLanz &: Hubeny J2003h f or the hot secondary. The 
model for the primary based on the parameters from lFremat et al. ( 2005 ) (see Table 3) makes 
an excellent match of the reconstructed primary spectrum. To find the best match of the re- 
constructed secondary spectrum, we created a grid of model UV spectra and cross-correlated 
each one with the reconstructed secondary spectrum. The largest amplitude ccf peak was 
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found at an interpolated effective temperature of 47.3 kK for models with log g = 4.00 and 
at 52.1 kK for models with log g = 4.75 (the highest value available in the published grid of 



Lanz fe Hubenyll2003l ). This (T e ff,logg) relation presumably makes the best overall match 
of the Fe IV to Fe V line depth ratios (and of those for other prominent ions). The ccf peak 
is slightly higher at the highest gravity, which is consistent with the appearance of the pres- 
sure broadened wings of He II A1640 (Fig. 5). We caution, however, that the models from 



Lanz fe Hubenyl (120031 ) assume solar abundances and a microturbulent velocity of 10 km s _1 , 
and changes about these assumptions will influence the derived best fit temperature. The 
model spectrum for T e g = 52.1 kK and log g = 4.75 is shown offset below the reconstructed 
secondary spectrum in Figures 3 to 7. 

Another comparison can be made to the observed UV spectrum of the sing l e, sdO star, 



BD+75°325. The spectrum of this star was analyzed in detail by lLanz et al.l (119971 ) who 
used the TLUSTY/SYNSPEC models to derive a helium abundance (He:H = 1 by number), 
T e ff = 52 kK, logg = 5.5, and a microturbulent velocity of 10 km s -1 . A comparison with 
BD+75°325 is particularly relevant because the hot companion in 59 Cyg is probably He- 
enriched due to past stripping by mass transfer. Fortunately, there are 92 UV spectra of 
BD+75°325 in the IUE archive that were made with the same SWP, high dispersion camera. 
We obtained these and processed them in the same way as done for the 59 Cyg spectra. 
Each spectrum was cross-correlated with a model TLUSTY/SYNSPEC spectrum to obtain 
the radial velocity (< V r >= — 59.4±0.3 km s _1 ), and then each spectrum was shifted to the 
lab frame before co-addition. Finally, the rectification of the mean spectrum was adjusted 
in the same way as done for the reconstructed secondary spectrum in order to inter-compare 
the spectra. The final UV spectrum of BD+75°325 is over-plotted in Figures 3 to 7, and, in 
general, it makes an excellent match to that of the hot secondary. 

The photospheric lines of both the hot secondary and BD+75°325 appear very sharp, 
and comparisons with model spectra indicate that any rotational line broadening is unre- 
solved (Usinz < 40 km s _1 ). Consequently, if we assume that the spectrum of BD+75°325 
is an effective match of that of the secondary, then we can find the monochromatic (mean 
UV) flux ratio /2//1 of the secondary to primary flux by adjusting the flux ratio used in the 
tomographic reconstruction so that the resulting line depths match those in the spectrum of 
BD+75°325. We derived this flux ratio by comparing the line depths in smoothed versions 
of the spectra in nine regions of ~ 50 A width that contained significantly deep lines. The 
best fit was obtained with a flux ratio f 2/ f total = 0.038 ± 0.005, and this value was used in 
the tomographic reconstruction shown in Figures 3 to 7. The formal error does not account 
for changes in line depth due to differences between secondary and BD+75°325 related to 
T e ff, logg, abundances, and microturbulence. Our impression from the good match made 
with the spectrum of BD+75°325 is that such differences are probably minor. Our results 
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are summarized in Table 3 that lists T e g, log g, and Vsini estimates for the hot subdwarf 
(column 3) along with those for the Be star (column 2) from lFremat et al.l (120051 ) . 



4. Heating of the Disk by the Companion 

One other means of detecting a hot companion is through obser vations of line emission 



from a heated region of the Be star disk that faces the companion (jHummel fc Stefll 12001 



20031 ). The localized heating effects are especially striking in the weak emission line of He I 
A6678 that usually displays a radial velocity curve similar to that of the companion but with 
a significantly larger semiamplitude. The He I A667 8 emission variation s are evi dent in the 
spectra of the known Be + sdO binari es, 59 Cyg (IMaintz et al.l 120051 ) , Per (jGies et al. 



19931 : IStefl et al.ll2000h . and FY CMa flPeters et al. 



20081). and their presence is the basis 



of th e proposed sdO companion of the Be star o Pup ( jKoubsky et al.l 12012c iRivinius et al. 
2012h . 



Here we revisit the He I A6678 emission line variations in the spectrum of 59 Cyg using a 
series of spectra we obtained with the Kitt Peak National Observatory C oude Feed Telescop e 
and additional spectra from the Database of Be Star Spectr E@(BeSS) dNeiner et alJboilh . 
Table 4 lists the origin, time span, spectral resolving power, number of observations, and 
observer name for the 59 spectra collected. All these spectra have a S/N ratio > 200, a 
resolving power greater than rs 10000, and record most of the He I A6678 feature (although 
in some of the BeSS spectra the blue or red wings are located at the boundary of the recorded 
spectrum). All the spectra were rectified to a unit continuum and were transformed to a 
uniform wavelength grid in the heliocentric frame. The spectra are presented in Figure 8 as 
a function of heliocentric velocity and orbital phase (from periastron at phase 0). The upper 
portion shows the individual spectral profiles with the continua aligned with the orbital 
phase, and the lower portion shows the spectral fluxes as a gray-scale image made through 
linear interpolation in orbital phase. The radial velocity excursions of the emission feature 
resemble those of the hot component (Fig. 1) but the emission has a larger semiamplitude 
(K em = 180 ± 10 km s" 1 ). The overal l appearance of the gray-scale figure is very similar to 
that presented by lMaintz et al.l (120051 ) (see their Fig. 2, which shows orbital phase increasing 
upwards), and this demonstrates that the emission line variations have continued through 
to the present time. 

The He I A6678 emission probably originates in the outer parts of the Be star disk that 
are directly illuminated by the flux of the hot subdwarf companion. The model presented 



2 http://basebe. obspm.fr/basebe/ 
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by iPeters et al.l ( 120081 ) for FY CMa (see their Fig. 8) can also be applied here to the case 
of 59 Cyg. In the reference frame of the Be star, the disk gas near the outer edge has a 
Keplerian velocity that exceeds that of the companion because of its location closer to the 
Be star. If we ignore third body effects, then we can estimate the approximate disk radius 
of line formation according to 

R em M x ( K 1 + K 2 \ 2 1 f K x + K 2 ~ 



a M 1 + M 2 V-fsTi + K cm J 1 + K x /K 2 \K X + K c 

where K\, K 2 , and K em are the semiamplitudes of the Be star, the subdwarf, and the He I 
emission, respectively. From the estimates given above, we obtain R era /a = 0.44, i.e., the 
emission region forms about midway between the Be star and subdwarf, presumably at a 
location near the outer boundary of the disk. 



5. Discussion 



The physical properties of the component stars are summarized in Table 3. The mass 
estimates are derived from t he orb ital results and an orbital inclination range of 60° to 80° 
advocated by iMaintz et al.l (120051 ) based upon the large projected rotational velocity and 
the past history of spectroscopic shell events. We can estimate the radii from the flux ratio, 
the spectral energy distribution, and the distance. The observed flux ratio is approximately 
related to radius ratio by 

h _ F2 ( R-i 
Ti~Ti\Ri, 

where F 2 /F% is the monochromatic flux ratio per unit area. Using the adopted temperatures 
in Table 3, we estimate that F 2 /F 1 = 10.6 at 1500 A based upon the TLUSTY/SYNSPEC 
models. Then the radius ratio is R 2 /Ri = 0.062 after correction for the small flux contribu- 



(2) 



tion of the Ab compon ent (assuming Am = —2.5 log /a&/(/i+/2) = 2.84; lMason et al.ll2009l ). 
Touhami et al.l (120121 ) estimated the angular diameter of primary Be star, 9 = 0.149 ± 0.004 
mas, by fitting the observed FUV spectrum to a composite model spectrum of Aa + Ab. 
This estimate should be considered as a value intermediate between the polar and equatorial 
diameters for the rapidly rotating Be star. Then to obtain the physical radius, we need to 
mul tiply the angular size times the distance. The original Hippa rcos distance was 34 5 ± 76 
pc (jPerryman et al.lll997l ). which was revised to 435 ± 79 pc by Ivan Leeuwenl (120071 ) . The 
radii given in Table 3 correspond to a distance range of 345 to 435 pc. 

Hot companion s of Be stars have now been d etected in three systems, 59 Cyg, FY CMa 
(IPeters et al.l 120081 ). and Per (IGies et al.l Il998l ). and there are some striking similarities 
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among these binaries. Table 5 lists the orbital period and eccentricity, mass, effective tem- 
perature, and luminosity for each component of these systems. All have a low mass ratio, 
M 2 /Mi « 0.11, and the temperatures of the respective components are very similar. Further- 
more, the hot sdO component is a slow rotator in each case, suggesting that the progenitor 
had realized internal synchronous rotation with the orbit prior to the completion of mass 
transfer. The hot subdwarfs in these binaries are more massive and brighter than the bet- 
ter known sd O stars in the field that represent a mixture of highly evolved, low mass stars 
(IHeberl 120091 ) . On the other hand, the subdwarfs in Be binaries have temperatures and lu- 
minosities that are si milar to those predicted by mass transfer models of binary evolution 
(jWellstein et al.ll200lf ). and their presence is strong evidence that these Be stars were spun 
up by an earlier stage of mass transfer when the stars were much closer. For example, if we 
assume conservative mass exchange, an in itial mass ratio of M 2 jM\ = 1.25 (typical of the 
systems modeled by IWellstein et al.ll200ll ). and a final mass ratio of M 2 /Mi = 0.11, then 
the current orbital periods should be ~ 21 x larger than their original periods. Thus, these 
systems probably began their lives as close pairs of B-stars. 

We also list in Table 5 the orbital eccentricities of the Be + sdO binaries. Only 59 Cyg 
has an eccentricity that is significantly different from zero, which is unusual because the 
intense tidal forces that the system experienced during large scale mass transfer would tend 
to circularize the orbit. We suspect that the massive third star in the system, Ab, may be 
a source of gravitational perturbations th at have led to a s ecular increase of the eccentricity 
following the conclusion of mass transfer (IFord et al.ll2000l ). 



The circumstellar disks of Be stars in Be + sdO systems are probably truncated by the 
tidal influence of the companion. We showed in Section 4 how the radius of the He I A6678 
emission region can be estimated based upon Keplerian motion (see eq. 1). The heated 
region probably forms near the outer disk boundary in the direction of the companion, and 
it is interesting to compare this emission radius with the Roche radius. We list in the final 
column of Table 5 the ratio of the emission radius R em ( iGies et al.l Il993l ; IStefl et al.l |2000| ; 



Peters et al.l 120081 ) to th e Roche radius at periastron (1 — e)i?R OC he (derived from the mass 
ratio and the formula of lEggletonl 119831 ) . We see that in all three cases the disk extends to 
at least 78 — 92% of the available Roche radius, so that in these active Be stars the disks are 
probably as large as permitted by tidal limitations. 

The example of 59 Cyg demonstrates that hot subdwarf companions of Be stars can 
be very difficult to detect. Indeed, it is thanks to the tenacity of the IUE observers that a 
sufficiently large number of UV spectra were available to make detection possible through 
the use of cross- correlation and Doppler tomography techniques. Thus, the question of how 
many Be stars may host hot subdwarf companions (and thus how many were spun up by 
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mass transfer) remains very uncertain because of the observational difficulties surrounding 
detection. Searching for such hot companions is of key importance to understand how Be 
stars became rapid rotators and to determine the numbers and kinds of their evolutionary 
progeny 
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was supported in part by NASA grant NNX10AD60G (GJP) and by the National Science 
Foundation under grant AST-1009080 (DRG). Institutional support has been provided from 
the GSU College of Arts and Sciences and the Research Program Enhancement fund of the 
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Table 1. IUE Radial Velocity Measurements 



Date UT SWP Orbital Vi ai (O - C)i V 2 cr 2 (O - C) 2 

(HJD-2,400,000) (yyyy-mm-dd) Number Phase (km s -1 ) (km s -1 ) (km s -1 ) (km s _1 ) (km s _1 ) (km s _1 ) 
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 



43708.914 1978-07-19 2050 0.156 -3.4 4.3 -15.1 

43864.202 1978-12-21 3664 0.665 12.9 3.2 17.7 59.9 6.7 -1.2 

44034.484 1979-06-09 5464 0.706 -11.3 3.8 -4.5 

44171.905 1979-10-25 6990 0.582 -13.3 3.2 -13.1 

44277.719 1980-02-08 7894 0.336 17.6 3.7 5.8 -118.6 5.6 -7.5 



Note. — A machine readable version of the full Table 1 is available in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 

I 

I— > 
Or 
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Table 2. Orbital Elements for 59 Cyg 



Element Harmanec et al. (2002) Maintz (2003) WE ccfs 



P (days) 28.1971 ±0.0038 28.192 ±0.004 28.1871 ± 0.0011 

T (HJD-2,400,000) . . • • • 51035.9 ± 0.2 45677.6 ± 0.3 

T sc (HJD-2,400,000) 50006.8 ± 0.6 51019.7 ±0.2 45693.0 ± 0.5 

ififkms- 1 ) 13.0 ± 1.0 24.8 ±3.9 11.7 ±0.9 

^(kms" 1 ) ••• 120.1 ±0.8 121.3 ±1.1 

7! (km s" 1 ) -16.5 ±0.8 ±4.3 ±1.5 ±2.1±0.7 a 

72 (km s' 1 ) ••• -0.9 ±0.5 -10.4 ±0.8 

e b 0.113 ±0.005 0.141 ±0.008 

wi (cleg) • • • 293 ± 3 257 ± 4 

M 2 /M l ••• 0.206 0.097 ±0.008 

M lS in 3 i(M ) ••• 7.23 6.08 ±0.14 

M 2 sm 3 i(M Q ) ••• 1.49 0.59 ± 0.05 

asini (R Q ) ••• 80.22 73.3 ± 0.8 

rmsi (km s _1 ) 3.1 ••• 8.1 

rms 2 (km s" 1 ) • • • • • • 9.1 



a Relative to the mean spectrum template. 
b Fixed. 



Table 3. Stellar Parameters for 59 Cyg 



Parameter 


Primary 


Secondary 


Teff(kK) 


21.8 ±0.7 


52.1 ±4.8 


log# (cgs) 


3.78 ±0.09 


5.0 ± 1.0 


V sin 2 (km s _1 ) 


379 ± 27 


< 40 


M/M & 


6.3-9.4 


0.62-0.91 


R/Rq 


5.5-7.0 


0.34 - 0.43 
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Table 4. He I A6678 Spectral Observations 



Source 


Dates 


R = A/AA 


Number 


Observer 


KPNO/CF 


1987 - 


1999 


25800 


22 


Peters 


BeSS/OHP/Elodie 


2003 


45000 


8 


Neiner 


KPNO/CF 


2004 - 


2008 


9500 


6 


Grundstrom 


BeSS/Castanet 


2008 - 


2012 


10000 


4 


Buil 


BeSS/Revel 


2009 


10000 


1 


Thizy 


BeSS/OVA 


2010 - 


2012 


10000 


18 


Mauclaire 



Table 5. Parameters for Be + sdO Binaries 



Star 


P 




log M/Mq 


log M/Mq 


logT cff 


logT efI 


log L/Lq 


log L/Lq 




Name 


(d) 


e 


(Be) 


(sdO) 


(Be) 


(sdO) 


(Be) 


(sdO) 


-Rem/[(1 — e)-RR,ochc] 


59 Cyg . 


28.2 


0.14 


0.9 


-0.1 


4.3 


4.7 


3.9 


3.0 


0.88 


FY CMa 


37.3 


0.00 


1.1 


0.1 


4.4 


4.7 


4.2 


3.3 


0.92 


Per ... 


126.7 


0.00 


1.0 


0.1 


4.4 


4.7 


4.3 


4.1 


0.78 
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Fig. 1. — The orbital phase variations of the cross-correlation functions with the hot star 
template spectrum after subtraction of the mean ccf. The ccf differences are shown as 
linear plots (top panel) and as a gray-scale image (lower panel). The intensity in the gray- 
scale image is scaled to ±4% and is calculated by a linear interpolation between the closest 
observed phases (shown by arrows along the right axis). The weak signal from the hot 
subdwarf spectrum appears as the "S" feature in the lower panel. 



- 19 - 




i i i i i i i i i i i 

0.0 0.5 1.0 

ORBITAL PHASE 



Fig. 2. — Radial velocity curves for the Be star and its companion. Orbital phase 0.0 
corresponds to periastron. The solid circles represent the ccf velocities for the Be star while 
the open circles show the same for the hot subdwarf. 
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Fig. 3. — A comparison of the O IV lines in the reconstructed UV spectrum of the secondary 
star (upper solid line) with those in the average spectrum of BD+75°325 (dotted line). Both 
are offset vertically by +0.5 to separate them from the model spectrum below (shown with 
no offset). All the spectra were smoothed by convolution with a Gaussian function with 
FWHM = 30 km s -1 (3 pixels). Line identifications for some of the stronger lines are given 
above the spectra. 
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Fig. 5. — A comparison of the He II line in the reconstructed UV spectrum of the secondary 
star {upper solid line) with those in the average spectrum of BD+75°325 {dotted line) in the 
same format as Fig. 3. The two lines in the lower part show how the He II line wings become 
narrower at logg = 4.00 (dashed line) compared to the nominal case with log g = 4.75 (solid 
line). 
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Fig. 6. — A comparison of the N IV and Fe IV lines in the reconstructed UV spectrum of the 
secondary star (upper solid line) with those in the average spectrum of BD+75°325 (dotted 
line) in the same format as Fig. 3. 
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Fig. 7. — A comparison of the N III lines in the reconstructed UV spectrum of the secondary 
star (upper solid line) with those in the average spectrum of BD+75°325 (dotted line) in the 
same format as Fig. 3. 
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Fig. 8. — The orbital phase variations in the He I A6678 emission line in the spectra of 59 Cyg 
are shown in linear plots (top panel) and as a gray-scale image (lower panel). The intensity 
in the gray-scale image is assigned according to the flux between the minimum (dark) and 
maximum (bright) observed values. The intensity between observed spectra is calculated by 
a linear interpolation between the closest observed phases (shown by arrows along the right 
axis). The scale bar at top right indicates the spectral flux relative to the local continuum 
flux. 



